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Abstract 


Eigenspace (ES) techniques are used to design an 
active flutter suppression system for the DAST ARW-2 flight 
test vehicle. The ES controller meets control surface 
activity specifications and at the flutter test condition 
provides reduced wing root torsion at the gust test condi- 
tion, and results in improved flutter boundaries. The 
ES controller is compared with a controller designed 
using Linear Quadratic (LQ) techniques. The LQ 
controller exhibits better phase margins at the flutter 
condition than does the ES controller but the LQ design 
requires large feedback gains on actuator states while 
the ES does not. This results in reduced overall actu- 
ator gain for the LQ design. 


Nomenclature 

Vectors 

u = control input 

v.= attainable closed loop eigenvector associated with X. eigen- 
1 value 1 

v d ^=desired closed loop eigenvector associated with eigenvalue 

w^= vector used in calculation of ES gain matrix, see Eqn. 17 

x = system state 


£ = flexure modal displacements 
£ c = control surface displacements 
T = disturbance input vector 


Matrices 


A 

B 



‘V 

[M s ] 


P. 

1 


[Q c ] 


R 

V 

w 

A 


aerodynamic coefficient matrix 
dynamics matrix 
control distribution matrix 
structural damping matrix 
aerodynamic coefficient error matrix 
control gain matrix 
structural stiffness matrix 
structural mass matrix 
eigenvector weighting matrix 
state weighting matrix 

calculated unsteady aerodynamic influence coefficient matrix 

s-plane representation of unsteady aerodynamic influence 

coefficient matrix 

control weighting matrix 

matrix whose columns are v. 

1 

matrix whose columns are w. 

1 

diagonal matrix composed of A ^ 


Scalars 


c = reference chord, 0.75m 
H(s)G(s) = loop transfer function 

= /IT 

= reduced frequency 

= reference length in gust model, 762m 
= Mach number 
= dynamic pressure 
= Laplace operator 
= forward velocity 
= aerodynamic lag frequencies 

= zero mean white noise input to gust model with intensity, 
I 


m 


k r 

V 


X. 

l 


= ith eigenvector 
a(E) = maximum singular value of E 

a(E) = minimum singular value of E 

oj = circular frequency 

€ g = normal'wind gust velocity 

= rms normal wind gust velocity 
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Subscripts 


i = inboard aileron 
o = outboard aileron 
ES = eigenspace 
LQ = linear quadratic 
s = structural 


Superscripts 

* = complex transpose 

T = transpose 
-1 = inverse 
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Eigenspace Techniques for Active Flutter Suppression 

I. Introduction 

The objective of the research described in this 
report is the application of Eigenspace (ES) design tech- 
niques to the synthesis of active flutter suppression 
systems. ES techniques allow the designer to use feed- 
back control to place closed loop eigenvalues and shape 
closed loop eigenvectors to satisfy performance specifi- 
cations. The basic theory behind ES design has been 
given by Moore [1] and others. Moore has shown that it 
is possible not only to place controllable eigenvalues 
and shape controllable eigenvectors but also to shape 
uncontrollable eigenvectors. Since the dynamic response 
characteristics of a system are determined by its eigen- 
vectors as well as its eigenvalues, the ability to shape 
eigenvectors provides the designer with an important 
tool. 

If performance specifications are given or can be 
interpreted in terms of desired closed loop eigenvalues 
and eigenvectors, ES techniques provide a natural design 
precedure where the desired eigenstructure (if obtainable) 
can be calculated directly without iteration. Cunningham 
[2] has shown how ES techniques can be used to improve 
aircraft flying qualities by placing rigid body poles 
in desired locations and decoupling rigid body modes. If 
performance specifications cannot be stated directly in 
terms of closed loop eigenvalues and eigenvectors, for 
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example specifications on rms responses or stability margins, 
the application of ES techniques is not so straightforward. 

This is the case in aeroelastic control problems such as 
flutter suppression. One of the principal contributions 
of this report is the description of a methodology for 
application of ES techniques to active flutter suppression 
and other aeroelastic control problems. 

To the authors' knowledge, the full power of ES tech- 
niques to shape eigenvectors as well as place eigenvalues 
has not been applied to active flutter suppression. Ostroff 
and Pines [3] used eigenvalue placement to design a flutter 
controller, but did not attempt to shape closed loop 
eigenvectors. In this report ES design techniques are 
applied to the design of a flutter suppression system for 
the DAST ARW-2 flight test vehicle. Only full-state feed- 
back is considered. Use of ES techniques in the design 
of robust state observers will be the subject of future 
investigations . 

It is shown that ES techniques can easily be used to 
design a flutter controller which satisfies performance 
specifications on rms control surface activity at the 
design condition. This controller also provides some 
gust load alleviation capability at off-nominal condi- 
tions. The ES controller requires no feedback on control 
surface actuator states and thus the open loop frequency 
response characteristics of the actuators are retained 
in the closed loop system. 
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The remainder of the report is divided into four 
major sections. First mathematical models of the air- 
craft, actuators, and normal wind gust are presented 
and performance requirements are discussed. Next the 
theory of ES design is described. Then the results 
obtained by applying this theory to the design of a 
flutter suppression system for the DAST ARW-2 aircraft 
are given and compared with results obtained from a 
controller designed using Linear Quadratic (LQ) optimal 
control theory. Finally conclusions and suggestions for 
future research are presented. A brief description and 
listing of a computer program used to perform ES design 
is given in Appendix B. 

II. Mathematical Models and Performance Requirements 

The DAST ARW-2 flight test vehicle is a Firebee II 
Drone which has been modified by replacing the conventional 
wing with a high aspect ratio, supercritical wing designed 
to flutter within the flight envelope. Two control 
surfaces, an inboard and an outboard aileron, are avail- 
able on the wing and a stabilizer is available on the hori- 
zontal tail. The outboard aileron is to be used for 
flutter suppression, gust load alleviation, and maneuver 
load alleviation. The inboard aileron is to be used for 
maneuver load alleviation, and the stabilizer is to be 
used to compensate for reduced static stability, for auto- 
matic flight control, gust load alleviation, and maneuver 
load alleviation. in practice only the outboard aileron 
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is to be used for flutter control. However, both the 
inboard and outboard surfaces were utilized in the controller 
designs since ES design methodologies are most useful in 
the design of controllers for systems with multiple control 
inputs . 

The design flight condition for the flutter control 
system was a Mach Number of 0.86 and an altitude of 4572 m 
(15000 ft). At this flight condition the flutter control 
system was required to stabilize the wing. a normal gust 
with rms velocity of 3.66 m/s (12 ft/s) was assumed at this 
condition. The rms deflection of the inboard aileron was 
limited to 10° and the deflection rate to 130 °/s. The rms 
deflection of the outboard aileron was limited to 15° and 
the deflection rate to 740°/sec. The control system was 
also to be evaluated on its ability to reduce bending and 
torsional stresses and shear forces in the wing at a gust 
test condition of Mach 0.7 and 4572 m. The vertical gust 
velocity at this condition was 18 m/s (59 ft/s) . In the 
actual DAST vehicle, the flutter suppression system would 
not be used at this condition because a separate gust load 
alleviation system is available. However, the gust test 
condition did provide a convenient condition for evaluating 
the performance of the flutter controllers at an off- 
nominal condition; therefore, in this study it was assumed 
that the flutter controller was also to be used for gust 
load alleviation. 
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Actuator/Control Surfaces 

Stabilizer . The stabilizer transfer function was a 
simple first-order lag 



The allowable rms control surface activity levels for the 
stabilizer were 7 degrees and 80 degrees/s. 

Inboard Aileron . An eleventh-order transfer function 
with third order numerator dynamics was given for the 
inboard aileron. Up to about 300 rad/s, a fourth order 
model gave an acceptable approximation of the eleventh 
order inboard aileron transfer function. This fourth order 
approximation was 

= 1.614X10 11 / (s 2 +67ls+(477) 2 ) (s 2 +322s+(878) 2 ) 

( 2 ) 

Outboard Aileron . A seventh-order transfer function 
with second-order numerator dynamics was given for the 
inboard aileron. Up to about 300 rad/s, a third order 
transfer function gave an acceptable approximation of the 
seventh order outboard aileron transfer function. This 
third order approximation was 

?C o/^ = 1.774xl0 7 / (s+180 ) (s 2 +251s+(314) 2 ) (3) 

Wind Gust . The wind gust was modeled by the second-oraer 
model given below 

S G /n(s) = (l+(/3L/V) s)/(l+(L/V) s) 2 (4) 
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Aircraft. The modes of the aircraft considered 


were rigid body plunge and pitch modes and seven symmetric 
aeroelastic modes of the wing. The flexible aircraft 
model was 


([M s ]s 2 + [C g ]s + [K s ]) U F ] + q[Q c (s) ] 


V 


= 0 (5) 


The aerodynamic influence matrix [Q c (s)] was calcul- 
ated over a range of reduced frequencies ranging from 
0.0 to 1.2 by a doublet lattice procedure. A rational s- 
plane approximation of Q c (s) was given by the approximation 


[Q a (s) ] 


[A o ] + [A 1 ] If + tV 


r cs 1 2 


2V 


L 
+ E 


tA m+2 ]s 


2v 

m=l s+— B 

c m 


This approximation has been widely used in design of active 
flutter suppression systems [4-10] . 

The first column of the matrix [A q ] was set equal to 
the first column of [Q (0)], which is zero since the aero- 
dynamic forces due to plunge displacement are zero. The 
second column of [A q ] gives the force due to a change in 
angle of attack resulting from a change in pitch angle and 

must be set equal to the second column of [Q (0) ] . The 

c 

first column of [A^] multiplied by c/2 gives the force 
due to a change in plunge velocity. Since the force due 
to a change in angle of attack must be the same regardless 
of whether this change results from a change in plunge 
velocity or a change in pitch angle, the first column of 
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[A.] must equal the second column of [Q (0)] divided by 
X c 

c/2 (there is also a scaling factor of 0.0062486 which 

must be divided into [Q (0)]). 

c 

Usual procedure is to select the B^s so as to 

bracket the reduced flutter frequencies [6]. Once the B's 

m 

are specified, the remaining elements of the [ A m ] matrices 
are determined to give the best least squares fit to Q c over 
the range of reduced frequencies for which this matrix has 
been calculated. 

In this study a new method was used to select the B's. 

An error matrix was defined as 

E ( j oj ) = [Q c ( jtu) ] - [Q a ( j oj) ] 

The norm of this matrix can be bounded above and below by 
its maximum and minimum singular values [11] 
g (E ) <ll E ^<g( E) 

The singular values of E are defined as the positive square 

roots of the eigenvalues of E*E. If the B's are chose to 

m 

minimize a(E), the norm of the error matrix will be small. In 

this study a single B was used in [Q ] • Since the reduced 

flutter frequency was known to be small (0.15), 8 was 

selected to minimize a( E) at zero frequency. The procedure 

for determining the [A ] ' s and 8 was as follows (1) a small 

m 

initial value for 8 was arbitrarily selected, (2) the first 

column of [A ] was set equal to zero, the second column of 

[A ] was set equal to the second column of [ Q c ( 0 ) ] and the 

first column of [A^] was set equal to the second column of 

[Q c (0)] divided by c/2 times a scale factor, (3) the remaining 

values of [A ] , [A.,], [A_] and [ A, ] were determined to give the 

O 1 2 3 
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best least squares fit to [Q c (j2kV/c)] over the range of 
reduced frequencies 0.0, 0.05, 0.2, 0.3, 0.3, 0.4, 0.5, 
0.6, 0.7, 0.8, 1.0 and 1.2 (4) the maximum singular 
value of the error matrix E(ju)=Q) was calculated (5) B 
was increased by a small. amount and the process was 
repeated until a(E) reached a minimum. The results are 
summarized below 


B a (E (0) ) 


0.100 

214.4 

0.125 

138.2 

0.130 

137.0 

0.135 

147.0 

0.150 

162.8 


Table 1 Singular Value of Error Matrix versus B 
The value of B equal to 0.13 which minimizes o’ ( E(0)) is 
very close to the reduced flutter frequency of 0.15. The 
eigenvalues resulting from this model for a Mach number 
of 0.86 and an altitude of 4572 m are given below 


Mode Eigenvalues 


plunge 

0, -1.17 

pitch 

-1 . 74± j6 . 29 

1st 

flexure 

41. 6±j 118.1 

3rd 

flexure 

-0 . 7± j 136 . 6 

4 th 

flexure 

-9 7 . 1± j 108 . 4 

6 th 

flexure 

-16 . 9 ± j 218 . 3 

8 th 

flexure 

-4 . 1± j 397 . 4 

9 th 

flexure 

-11 . 0 ± j 4 25 . 0 

10th 

flexure 

-24 . 8± j 452 . 5 


Table 2 Eigenvalues of DAST ARW-2 Flight Test 
Vehicle at M=0.86, h=4572:n based on 2-Rigid Body 
and 7 Flexure Modes 
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The values of the flexure mode eigenvalues are close to 

those calculated for the ARW-2 using a model containing 

four lag states [12] . The pitch eigenvalues were not 

too different from the actual values of -1.5±j5.6; 

however, one of the plunge eigenvalues is considerably 

different from the actual values of 0 . 0 ± j 1 . lxlO -3 . The 

flutter characteristics of the aircraft are not affected 

by the rigid body modes and these modes are not included 

in the model used in the control system design. The 

poorest correlation between elements of the [Q ] and [ Q ] 

C A 

matrices occurred for the coefficients associated with the 

gust velocity, even then rms structural and control surface 

gust responses calculated using the model with a single B 

were close to those calculated from higher-order models 

using several values of B * 

m 

It is felt that selection of the numerical values of 

the B's based on the maximum and minimum singular values of 
m 

an error matrix such as given by Eq. 7 has considerable 
potential in generating low-order approximate models of 
unsteady aerodynamics. Even the simple approach of mini- 
mizing the maximum singular value of the error matrix at 
a fixed frequency yielded acceptable results. Other 
approaches such as choosing the B^s to minimize functions 
of both the maximum and minimum singular values over a 
range of frequencies should be investigated. 
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Fig. 1 Root Locus of Seven Mode Model of Uncontrolled Wing 
as Velocity is Varied 
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Some possible indices are as follows: 

1. Minimize average a and £ over a range of frequen- 
cies, i.e. 

, 0) 2 

min{- — — - / c (E) ( joi) ) do)} 

0 ), 

2 1 o) j 

2. Minimize maximum value of a and £ over a range of 
frequencies, i.e. 

min {max (a(E(jo)) ) } o) x £ o) £ o) 2 

3. Minimize weighted sum of a and £ over a range of 

frequencies, i.e. 

N 

min{ L a . a (E ( jo) . ) ) 
i=l 1 

An examination of the 10 flexure modes (see Appendix A) 
indicated that modes 2 and 5 were primarily fuselage bending 
modes and mode 7 was exculsively a tail mode. Therefore 
these three modes were not considered further in the anal- 
ysis. Mode 1 was primarily the first wing bending mode, 
mode 2 was the second wing bending mode, and mode 6 was 
the first wing torsion mode. These modes were obviously 
important in modeling flutter and were retained. Modes 
3 and 8 included wing tip bending and it was felt that 
these modes should also be analyzed further. Mode 9 was 
primarily wing bending and mode 10 was primarily wing torsion 
and these modes were also retained. Thus seven flexure 
modes, the 1st, 3rd, 4th, 6th, 8th, 9th and 10th were 
included in the initial flutter analysis. The loci of the 
eigenvalues of these modes as velocity was varied are shown 
in Fig. 1. The first mode flutters at a velocity of about 
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241 m/s (M=0.75) at a frequency of 120 rad/s. Modes 3 and 

8 are insensitive to changes in velocity indicating 

that they are primarily vibrational modes; however, modes 

9 and 10 do vary somewhat with velocity. Eigenvalues 

were calculated using models in which various modes were 

deleted. Deletion of the 3rd, 8th, 9th and 10th modes 

has very little affect on the lower modes. 

3 Mode Model 

Mode No Eigenvalue 

6 -15 . 7± j 216 . 1 

4 -94 . 3± jl06 . 4 

1 +39 . 9± j 118 . 0 

Table 3 Eigenvalues of DAST ARW-2 Flight Test Vehicle 
at M=0 . 86 h=4572m Based on 3 Flexure Modes. 

Basic flutter characteristics could be accurately 
modeled with only three modes corresponding to the first 
and second bending (modes 1 and 4) and first torsion 
(mode 6) , and computational expense could be reduced signif 
icantly by reducing the order of the system model* there- 
fore, the design studies were based on a structural model 
containing only these three flexure modes. The loci of 
the eigenvalues of the three mode model as velocity is 
varied are shown in Fig. 2. 

Equations 2-6 can be combined to give the equations 
which describe the wing, control surfaces and actuators, 
and wind gust in vector-matrix form as 

x = Ax + Bu + Tn (8) 

The 18th order state vector, x, consists of (1) the dis- 
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H 


Fig. 2 Root Locus of Three Mode Model of Uncontrolled Wing 
as Velocity is Varied , 
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placements and velocities associated with the three flexure 
modes, (2) the unsteady aerodynamic lag states, (3) the 
states associated with the inboard and outboard control 
surfaces, and~(4) the states associated with the wind gust 
model. The 2nd order control vector, u, consists of the 
inputs to the inboard and outboard aileron/actuators and 
the zero mean white noise input, n, drives the wind gust 
model. (See Appendix B for a more detailed discussion of Eq. 8) 

III. Eigenspace Design 

Theory 

Moore [1] and others have shown how feedback can be 
used to directly place eigenvalues and also shape eigen- 
vectors. If performance specifications are given or can be 
interpreted in terms of desired closed loop eigenvalues 
and eigenvectors, then ES techniques can provide a natural 
design procedure. If performance specifications cannot be 
clearly stated in terms of closed loop eigenvalues and 
eigenvectors, for example specifications on rms responses, 
it may be necessary to iterate eigenvalues and eigenvectors 
until performance specifications are met. 

Before discussing the details of the ES flutter con- 
troller, the. theoretical basis of ES design methodology 
will be discussed. Consider a system modeled in state 
variable form as 

x = Ax + Bu + Tn 
with measurements 
y = Cx 
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where the state x£R n , the feedback control u£R m , and the 
output y£R°. If we express u in output feedback form 
u = Ky = KCx 

the the closed loop system response is 

x = (A + BKC) x + Tn (9) 

The ES design procedure consists, of determining a 
gain matrix, K, such that for all desired closed loop 
eigenvalue (A^) and eigenvector (v^) pairs 

(A + BKC) v . = A . v . (10) 

l li 

Introduce wtR where 

w. = KCv i (11) 

The problem becomes one of determing w^ such that 


(A. I-A) v. = Bw. 
ill 

where K is determined from Eq. 11. 


Case 1: Full State Feedback 

If all the states of the system are available, the 
feedback control law becomes (assuming p=n and letting 
C=I without loss in generality) 
u = Kx 

The design problem is to find w.^ and K such that 

(A.I-A)v. = Bw. M2) 

ill ' ' 

and 


w. = Kv. 

l l 

For the entire collection of closed 
eigenvectors Eq. (12) becomes 
VA - AV = BW 


(13) 

loop eigenvalues and 

(14) 
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Likewise from Eq. (13) 

W = KV (15) 

therefore 

K = WV -1 (16) 

Case 2: Output Feedback (p<n) 

y = Cx 

In this case Eq. (15) becomes 
W = KCV 

Since C is of rank p<n a unique solution for K is impossible. 
The alternative is to place only p eigenvalue-eigenvector 
pairs, i.e. choose p X^s and the corresponding p columns of 
V and solve for K such that 
K = W(CV ) -1 

The key design issue is to find W which satisfies Eq. (14) . 

In general one cannot completely satisfy both exact eigen- 
value and eigenvector placement. 

Case 3: Single Input Systems 

For single inputs Eq. (12) reduces to 
( X . I- A) v . = bw . 

l li 

where w^ is a scalar. This single variable cannot be 
adjusted to place n parameters on the left hand side, there- 
fore, the single input case only involves pole placement 
with arbitrary eigenvector position. The gain K is unique 
in this case. 

Case 4: Multiple Input Systems 

The real benefit of eigenspace techniques is realized 
when more than one control is available. If the rank of B 
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is n, i.e., one independent control for each state then 
from Eq. (12) 

w. = B -1 ( A . I -A) v. 
i 11 

for each desired A^ and v^ of the closed loop system. 


Practically speaking, however, one has fewer controls than 
states and exact placement of A^ and for all i=l,2,...n is 

impossible. The design procedure then becomes one of choosing 
portions of v_^ to eliminate certain state responses from a 
mode while emphasizing others and letting other responses 
react arbitrarily. Assuming the system is controllable and 
rank (B) = m<n, only m free parameters can be specified. If 
it is desired to change an eigenvalue associated with a control- 
lable mode, then for the new eigenvalue, A^ the matrix 
(A^I-A) is nonsingular. Thus 

v. = (A.I-A) _1 Bw. (17) 

li i 

or 

v . = L . w . (18) 

l li 


where 


L i = (X i I_A) ~ 1b 


Since there are not enough independent controls available 

to arbitrarily place all A^'s and v^'s, the w^'s are selected 

to minimize the following least square performance index 

J. = (v, -v. )*P.(v. -v.) 
i d . l i d . i 

l i 

where P^ is used to emphasize certain components of v^ , . 
Solving for w^ which minimizes produces an optimal psuedo- 
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inverse 


w. = (L*.P.L.) -1 L*.P.v, 

1 111 lid. 

1 

and is given by Eq. (18) . 

If certain eigensolutions are not to be altered (e.g. 
actuator poles) or if X. is uncontrollable then 
w . =0 

l 

and Eq. (12) is satisfied. The general design procedure 

used is to vary P., X., and v d until performance specif i- 

1 i 

cations on rms values of the state and control are met. 

The above analysis is valid for complex or real eigen- 
solutions. If an entirely real K matrix is desired then 
the the problem must be decoupled into its real and imaginary 
parts. If 

1 = 1 R + jlj 

V = V R + jVj 
w = w R + jw I 
then for real A 


(X.I-A)v. = Bw. 
i i l 

becomes in entirely real terms 


( XrI-A) -XjI 

V R 


l 

s? 

1 

Xji ( XrI-a) 

_ V I_ 


_ BW ! J 


(19) 


and 


V R 


(XrI-A) -X i I 

J 1 - 


xi (XrI-a) 



( 20 ) 
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or 


where 


U r I-A) -XjI 

XjI (A R I-A) 




Then the optimal psuedo-inverse solution in entirely real 


terms is 


where 


T -IT 

w. = (L. p.L.) L. p . V 
1 111 11 d, 


W i = W R . 


v d. = v d R 

l R i 




P . = 

l 


P R. 

l 


To determine K with A^ and A 2 complex conjugates. 


we must solve 


K[v r +jv ,v -jv ,V] - [w +jw ,w -jw ,W] 

1 1 1 1 R 1 T 1 R 1 h 

where V and W are the remaining ' s and w^'s. Post 
multiplying both sides of this equation by 


h -jh 
R = h +jh 
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yields 

K[v p #v,,V] = [w R »w T »W] (25) 

R 1 I 1 R 1 Z 1 

Now the left hand eigenvector matrix is nonsingular so 

K = ,w T , W] [v , v , V] 1 (26) 

R 1 I 1 R 1 I 1 

This procedure can be applied to any complex conjugate 
pair. It should be noted, however, that the transforma- 
tion matrix R is not unique so neither is K. 

Application to Active Flutter Suppression 

Using the above procedure a flutter control system 
was designed. The initial ES controller was designed 
by rotating the unstable eigenvalues around the imaginary 
axis and leaving all other eigensolutions in their open 
loop configuration. This resulted in acceptable rms control 
surface activity at the flutter condition and a stable 
response at the gust test condition. Although the ES 
design stabilized the wing at the gust test condition, the 
maximum allowable values for the rms inboard deflection 
rate and outboard deflection and deflection rate were exceeded. 
It was felt that the performance at the gust test case 
might be enhanced by redesign of the control system. Since 
the aircraft exhibits satisfactory response at velocities 
somewhat less than the flutter speed, it was decided to use 
ES design techniques to force the closed loop eigenvalues 
at the design velocity (275 m/s) to be the same as the 
open loop eigenvalues at the velocity of 200 m/s (this is 
20% less than the flutter speed) and to force the closed 
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loop eigenvectors at 275 m/s to approach the open loop 
eigenvectors at 200 m/s. The control was designed to retain 
the eigenvalues and eigenvectors associated with the actu- 
ators, unsteady aerodynamics, and gust model at their open 
loop values so that energy would not be transferred to these 
modes. (The gust states are uncontrollable and the gust 
eigenvalues cannot, in any case, be moved.) The ability to 
keep poles in desired locations is one of the key advantages 
of the ES design. 

In the open loop condition, the actuator dynamics are 
decoupled from the aeroelastic modes; therefore, the compon- 
ents of the aeroelastic eigenvectors in the direction of the 
actuator states are zero. Thus driving the closed loop 
aeroelastic eigenvectors to exactly their open loop values 
would decouple the actuators from the aeroelastic response 
of the wing. This is obviously not desirable since the 
wing would be uncontrolled; however, it is possible to 
reduce control surface activity and still stabilize the wing 
by using the weighting matrices to penalize large values 
of the components of the closed loop aeroelastic eigen- 
vectors in the actuator directions. Initially all closed 
loop aeroelastic eigenvalues at 275 m/s were placed at the 
locations of the open loop aeroelastic eigenvalues at 200 m/s 
and the weighting matrices were chosen to be identity 
matrices. This resulted in a reduction of control surface 
activity of less than 7% at the gust test condition. The 
rms inboard aileron deflection rate was still over three 
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times its allowable value while the outboard aileron deflec- 
tion and rates were about twice their allowable values. A 
gust of 5.7 m/s would saturate the inboard control surface 
while the outboard surface remained unsaturated since its 
deflection and rate were about two-thirds of the allow- 
able values. 

It was decided to shift some of the control effort from 
the inboard to the outboard aileron in an attempt to increase 
the gust velocity at which the system would saturate. This 
was accomplished by increasing the weights on the components 
of the aeroelastic eigenvectors in the inboard actuator 
directions while retaining all other weights at one. Values 
of these weights were increased to 2.5xl0 3 yielding the results 
in Table 4. The inboard rate was reduced substantially 
without excessively increasing outboard activity. Specif- 
ications on inboard rate and outboard deflection and rate 
were still not met at the gust test condition but all of 
the quantities were about twice their maximum values . Thus 
both control surfaces would saturate at about the same gust 
velocity (8.3 m/s). Since it proved impossible to further 
reduce outboard surface activity by adjusting the weighting 
matrices P^, the ES controller resulting from eigenvector 
shaping to reduce inboard rate was chosen for further eval- 
uation. 
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Table 4 


Comparison of RMS Control Surface Activity 
for Various Controller Designs and Flight Conditions 



Inbd Defl 

WtiSZ&M 

Inbd Rate 

(Dea/s) 

Outbd Defl 

(Deg) 

Outbd Rate (Deg/s) 

Controller 

Design 

Flutter 

Gust 

Test 

Flutter 

Gust 

Test 

Flutter 

Gust 

Test 

r.7 J.J. Gust 

Flutter Test 


Unstable Roots 
Rotated About 
Imag. Axis 


ES 

0.5 

5.0 

108.0 

412.0 

3.3 

31.4 

509.0 

1464.0 

Unstable Roots 
Rotated About 
Imag. Axis 

LQ 

1.8 

* 

271.0 

* 

2.8 

* 

407.0 

* 

Eigenvector 
shaping to 
Reduce Inbd 
Rate (Final ES) 

0.9 

8.1 

86.0 

279.7 

4.7 

31.2 

612.0 

1464.0 

Weighting 
Inbd Rate 
in Perf Index 
(Final LQ) 

0.7 

4.6 

85.0 

284.1 

3.1 

30.2 

486.0 

1335.1 


Final ES 
Inbd Failed 

- 

- 

- 

- 

4.4 

26.5 

572.0 

1427.0 

Final LQ 
Inbd Failed 

- 

- 

- 

- 

3.5 

30.0 

519.0 

1410.0 


Max Allowable 

10.0 

10.0 

130.0 

130.0 

15.0 

15.0 

740.0 

740.0 


*Uns table 
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Results 


In addition to the ES design, a flutter controller 
was designed using Linear Quadratic (LQ) optimal control 
theory. LQ control theory 'has been discussed extensively 
in numerous texts, for example Ref. 13. The basic LQ 
design procedure consists of selecting quadratic weighting 
matrices Q and R in a scalar performance index 

J = H /°°[x T Qx + u T Ru] dt (2 7) 

o 

The control which minimizes this performance index is given 
as 

u = Kx 

where K is the solution of a matrix Ricatti equation. 

It is well known that if Q is set equal to zero and R 
equal to the identity matrix in the performance index given 
by Eq. (27), then all stable eigenvalues will remain 
unchanged and all unstable eigenvalues will be rotated 
about the imaginary axis [13]. Initial LQ design was 
performed using this approach; however, as shown in Table 
4, the rms deflection rate for the inboard actuator was 
approximately double its allowed maximum at the flutter 
test condition and the controller resulted in an unstable 
wing at the gust test condition (note the uncontrolled 
wing is stable at this condition) . Since all rms responses 
except the inboard actuator deflection rate were acceptable 
at the flutter test condition, only this state was weighted 
in the performance index. The results of varying the 
weight on inboard actuator deflection rate is shown in Table 
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Table 5 


RMS Response for Various Weights 
on Inboard Actuator Rate 


Weight on Inbd 
Deflection Rate 

RMS Response (Degrees and Degrees/s) 

Inbd Def. 

Inbd Rate 

Outbd Def. 

Outbd Rate 

0 

1.60 

217.6 

2.73 

427 

1 x 10 -5 

1.58 

256.0 

2.74 

428 

5 x 10~ s 

0.92 

122.0 

3.02 

470 

1 x 10“ 4 

0.66 

84.8 

3.14 

486 

5 x lO -4 

0.22 

26.0 

3.32 

514 

1 x 10~ 3 

0.12 

14.0 

3.36 

518 

Max Allowable 

10.0 

130.0 

15.0 

740 
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5. The rms value of the inboard actuator deflection rate 
decreases fairly rapidly as its weight is increased and the 
rms activity level of the outboard control surface does 
not increase substantially. The eigenvalues associated 
with the flexure modes and the outboard aileron do not 
change as the weight on the inboard actuator rate is 
changed; however, as shown in Fig. 3, the eigenvalues asso- 
ciated with the inboard actuator change radically. The 
moduli of three of the roots become very large while the 
fourth root approaches zero. Since any large change in 
actuator frequency response characteristics would be diffi- 
cult to obtain without substantially redesigning the actu- 
ator, a value for the weight on actuator deflection rate 
-4 

of 1 x 10 was selected. This gave acceptable rms responses 
at the flutter test condition and also stabilized the wing 
at the gust test condition and did not affect the frequency 
response of the inboard actuator too much. The open and 
closed loop inboard actuator frequency response curves for 
the final design are shown in Fig. 4. It can be seen that 
the overall gain for the closed loop system is reduced 
compared with the open loop response. At zero frequency 
the open loop gain is unity and the closed loop gain is 0.79. 
Thus the forward loop gain would have to be increased by 

21% in order to restore the steady state gain to its open 
loop value. This might be difficult without modifications 
to existing actuator hardware. 
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Fig. 3 Variation in Actuator Root 
with Weighting on Inboard Aileron 
Deflection Rate 
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Fig. 4 Inboard Actuator Transfer Function Frequency Response 


It appears to be impossible to design either an ES or 
LQ flutter controller which also meets the specification on 
rms control surface activity at the gust test condition. 

In the actual DAST ARW-2 vehicle, the inboard aileron is 
used for manuever load alleviation but not for flutter sup- 
pression or gust load alleviation. Furthermore a different 
controller design is used for gust load alleviation than for 
flutter suppression. Thus it is not too surprising that the 
ES and LQ flutter controllers do not meet specifications on 
rms surface activity at the gust test condition. However, 
if the rms gust velocity at the gust condition is reduced to 
8.3 m/s (about 50% of its nominal value) both the ES and LQ 
controllers meet all specifications on rms surface activity. 
Table 6 shows that both controllers provide some torsional 
load alleviation at this flight condition. The bending and 
torsional stresses and shear force were calculated at nine 
stations on the wing (station 1 is at the root and station 
9 is near the tip) . The LQ design results in the largest 
rms bending moments at all stations. The ES design results 
in bending moments which are lower than those resulting 
from the LQ design but higher than the uncontrolled values. 
The differences in bending moment in all three cases are 
not large. Both the LQ and ES designs result in substantial 
reductions (over 50%) in torsional stresses at all stations 
compared with the uncontrolled values. The LQ reduces 
torsional stresses slightly more than the ES design. The 
ES design results in lowest values of shear near the wing 
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Table 6 


RMS Loads at Modified Gust Test Condition 
(M= . 7 , h=4572m, V G =8 . 3 m/s) 




Root 



Station 




Tip 



1 

2 

3 

4 

5 

6 

7 

8 

9 


Bending 

Moment 

(N-M) 

5740 

4530 

3426 

2444 

1371 

889 

377 

83 

5 

FINAL 

LQ 

Shear 

(N) 

92 

82 

74 

63 

42 

36 

24 

6 

1 


Torque 

(N-M) 

36 

54 

73 

82 

73 

63 

52 

6 

2 


Bending 

Moment 

(N-M) 

5558 

4412 

3359 

2418 

1374 

900 

389 

91 

5 

FINAL 

ES 

Shear 

(N) 

87 

79 

71 

61 

42 

36 

24 

6 

1 


Torque 

(N-M) 

46 

63 

79 

85 

76 

65 

53 

7 

1 


Bending 

Moment 

(N-M) 

5545 

4262 

3116 

2125 

1115 

671 

248 

47 

3 

OPEN 

LOOP 

Shear 

(N) 

98 

85 

75 

61 

38 

31 

18 

3 

1 


Torque 

(N-M) 

106 

130 

152 

156 

139 

121 

101 

24 

3 
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root but the uncontrolled values of shear are less near 
the tip. The differences in shear between the LQ, ES, and 
uncontrolled cases are not large. 

The root locus of the eigenvalues of the wing with the 
LQ controller is shown in Fig. 5 as velocity is varied. 

The wing goes unstable at a velocity of about 295 m/s (the 
design condition was 275 m/s and the uncontrolled flutter 
speed was 241 m/s) . It is interesting to note that for the 
LQ controller, one of the roots associated with the unsteady 
aerodynamic lag states goes unstable resulting in zero 
frequency flutter. In the uncontrolled case, the first 
bending model goes unstable in classical coupled mode flutter 
(Fig. 1). The root locus for the ES design is shown in Fig. 6. 
For the ES controller the wing goes unstable at a velo- 
city of about 310 m/s. As in the uncontrolled case the 
first bending mode goes unstable, but in the controlled case 
the eigenvalues associated with this mode move to the real 
axis where one real root goes unstable resulting in zero 
frequency flutter. 

The results of varying altitude while maintaining Mach 
number constant are shown in Fig. 7. At M=0.86, the uncon- 
trolled wing is unstable until an altitude of 6700 m is 
reached. At the same Mach number, the LQ controller results 
in a stable wing at all altitudes above 3200 m and the ES 
controller stabilizes the wing above altitudes of 2900 m. 

At M=0 . 7 the uncontrolled wing is stable for altitudes 
above 1800 m whereas the ES controller stabilizes the 
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Fig. 6 Root Locus of ES Controlled Wing as Velocity is Varied 
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Fig. 7 Flutter Boundaries for Uncontrolled, LQ, and ES 
Controlled Wing 
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wing for altitudes above 2100 m and the LQ controller 
stabilizes the wing for altitudes above 2400 m. 

The ability of the ES and LQ systems to stabilize the 
wing in case of an actuator failure yielded little difference 
in the performance of the two systems. In the case where 
an inboard actuator failure was simulated, the outboard 
aileron was capable of stabilizing the wing with only small 
increases in rms surface activity (see Table 4) . When 
the outboard actuator failed the inboard aileron was unable 
to stabilize the wing in both designs. 

Since the outboard aileron is critical in stabilizing 
the wing, it is important to examine the stability margins 
associated with the outboard control loop. This can be 
accomplished by examining the loop transfer function 
H ( s) G ( s) = K T [Is-A] -1 B 

Note since u = +Kx, the characteristic equation is given 
by 

1 - G(s) H (s) = 0 

and the critical condition occurs when the phase angle of 
H ( G ( jw) is zero. For the LQ design 


H (s) G (s) 


-148.8 (s-2.6) 


LQ (s z -79.8s+(124.5) 2 ) 

It is interesting to note that the actuator poles, the 
poles associated with the unsteady aerodynamic lag states 
and all of the poles associated with the stable flexure 
modes have been cancelled by zeros. The only remaining 
poles are those associated with the unstable first bending 
mode. Bode diagrams for the LQ loop transfer function are 
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given in Figs. 8 and 9 . The gain margin is 6Db and the 
phase margins are 60°. This is not surprising, however, 
since Safanov and Athans [14] show that LQ controllers 
yield excellent gain and phase margins. 

The loop transfer function for the ES controller is 


H (s) G (s) Es 


5 1 . 4 5 ( s-4 2 7 . 4 ) (s-66.5) (s+4 2. 4 ) 
(s‘-79.8s+<124.5)‘) ( s ‘+18 8 . 6 s+ ( 142 . 2 ) * ) 


Again all actuator poles and poles associated with unsteady 
aerodynamic lag states have been cancelled by zeros. But 
only the poles associated with the first torsion mode 
have been cancelled and the unstable first bending mode and 
stable second bending mode remain. Bode diagrams for this 
transfer function are also shown in Figs. 8 and 9 . Gain 
margins are 6 Db or better but phase margins are less than 
20°, Gain versus frequency plots for the LQ and ES transfer 
functions are very similar for frequencies above 10 rad/s. 
Both curves peak at the flutter frequency with a gain of 
6 Db and then roll off with increasing frequency. 


Conclusions 

The ES and LQ controllers give very similar results in 
terras of required control surface activity. At the gust 
test condition the ES design exhibits lower wing root bending 
moment and shear than the LQ design but the LQ controller 
provides slightly lower torsional stress. Both the ES 
and LQ designs provide significantly reduced torsional 
stress at the wing root, slightly reduced shear, and' 
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slightly increased bending moment compared with the open 
loop response. Both the LQ and ES controllers significantly 
increased flutter speed compared with the uncontrolled 
wing. The ES controller results in a slightly greater 
flutter speed than the LQ controller. For a fixed Mach 
number the ES design is stable at lower altitudes than the 
LQ design. The LQ design exhibits significantly better 
phase margins than the ES design at the flutter test condi- 
tion. The gain margins are the same. Since the phase 
margins are determined near the flutter frequency, it is 
very important that the ES design be based on a model which 
is accurate in this frequency range. Both the LQ and ES 
designs exhibit excellent roll off at higher frequencies 
where modeling uncertainties are large. 

The LQ design requires large feedback gains on the 
inboard actuator states. This reduces the overall inboard 
actuator gain. Increased forward loop gain would be required 
to restore open loop characteristics. This might necessi- 
tate redesign of existing actuator hardware. The ES controller 
does not require actuator feedback and the closed loop 
frequency response characteristics of the actuators is the 
same as the open loop case. Since the inboard actuator is 
not very effective for flutter control this is not an impor- 
tant consideration for the ARW-2, but could be critical in 
other applications. The ES controller requires only the 

structural and aerodynamic states, thus a lower-order 
observer could be used to realize this controller than would 
be required by the LQ controller. Finally the computational 
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algorithms required to calculate the ES controller gains 
are simpler and less expensive to use than those required 
to calculate the LQ controller gains. 

The authors are currently working on ES design techni- 
ques that shape eigenvectors associated with uncontrollable 
states. It appears to be necessary to shape the eigen- 
vectors associated with the gust states in order to further 
reduce rms control activity. Work is also in progress on 
improving ES stability margins and realization of ES control- 
lers by means of observers and direct output feedback. 

Doyle and Stem [15] have shown that if a Kalman Filter 
is used for state estimation, the robustness properties of 
full-state feedback can be recovered by introducing ficti- 
tious plant noise in the Ricatti equation used to obtain 
the filter gain matrix. As the magnitude of this plant 
noise approaches infinity, the filter poles asymptotically 
approach the transfer zeros (if the transfer zeros are in 
the right half of the complex plane the filter poles approach 
the mirror image of the zeros in the left half plane) or 
approach infinity in a Butterworth pattern. Since the Doyle- 
Stein procedure gives the location of the filter poles, it 
is possible to directly obtain the gain matrix which yields 

these poles using ES design techniques without solving the 
filter Ricatti equation. In the single-output, single- 
input case, this should yield the same result as obtained 
from solving the Ricatti equation; however, in the multi- 
input, multi-output case there is not a unique gain matrix 
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which yields a specified set of poles; thus, ES techniques 
will necessarily not yield the same results as obtained 
by the Doyle-Stein procedure. The utility of ES techniques 
in observer design should be studied in more detail. In 
addition, the effects on system performance of including 
flexure modes that were neglected during the design phase 
needs to be studied. 


Acknowledgement 

The results on ES design were obtained under Grant 
NAG-1-217 from NASA Langley Research Center. Mr. William 
Adams of NASA Langley and Mr. San jay Garg of the University 
of Minnesota were instrumental in developing and interpreting 
the mathematical models of the ARW-2. Mr. X.G. Li of the 
University of Minnesota played an important role in calcu- 
lating the results for the LQ control system. Thanks are 
also due Dr. Thomas B. Cunningham of Honeywell, Inc. for a 
number of useful discussions concerning Eigenspace tech- 
niques . 


- 46 - 



References 


1. Moore, B.C., "On the Flexibility Offered by Full- 
State Feedback in Multivariable Systems Beyond Closed 
Loop Eigenvalue Assignment," IEEE Trans. Auto. Control , 

Vol. 21, pp 682-691; 1976. 

2. Cunningham, Thomas B., " Eigenspace Selection Procedures 
for Closed Loop Response Shaping with Modal Control ," 
Proceedings IEEE Conference on Decision and Control; 

Dec. 1980. 

3. Ostroff, A . J . and Pines, S., "Application of Modal 
Control to Wing - Flutter Suppression ,” NASA Tech 
Paper 1983; May 1982. 

4. Newsom, J.R., "Control Law Synthesis for Active Flutter 
Suppression Using Optimal Control Theory," Journal of 
Guidance and Control , Vol. 2; Sept. - Oct. 1979 pp 388-394. 

5. Abel, I., Newsom, J.R., and Dunn, H.J., "Application 

of Two Synthesis Methods for Active Flutter Suppression 
Using Linear Quadratic Gaussian Control Theory ," AIAA 
Paper 79 - 1633; Aug. 1979. 

6. Mahesh, J.K. , Stone, C.R., Garrard, W.L., and Dunn, H.J., 

"Control Law Synthesis for Flutter Suppression Using 
Linear Quadratic Gaussian Control Theory," Journal of 
Guidance and Control , July - Aug., 1981, pp 415-422. 

7. Mukhopadhyay , V., Newsom, J.R., and Abel, I., "Reduced- 
Order Optimal Feedback Control Law Synthesis for Flutter 
Suppression," Journal of Guidance and Control , July - 
Aug. 1981, pp 382-395. 


- 47 - 



8. Gangsaas, D., Ly, U. , and Norman, D.C., "Practical 
Gust Load Alleviation and Flutter Suppression Control 
Law Based on LQG Methodology , " AIAA Paper 81 - 0021, 

Jan. 1981. 

9. Gangsaas, D., and Ly, U. , "Application of a Modified 
Linear Quadratic Gaussian Design to Active Control of 

a Transport Airplane," AIAA Paper 79 - 1746, Aug. 1979. 

10. Dowell, Earl H. , "A Simple Method for Converting Fre- 
quency Domain Aerodynamics to the Time Domain," NASA 
Tech Memo 81844, Oct. 1980. 

11. Noble, Ben, and Daniel, James W., "Applied Linear 
Algebra," 2nd Edition, Sections 9.6 and 11.6, Prentice 
Hall, 1977. 

12. Adams, William M. , Jr. and Tiffany, Sherwood H., "An 
Updated Comparison of NASA/Boeing Predictions of the 
Open Loop , Symmetric Flutter Characteristics of the 
DAST AEW-2," APCO TM 81-1, NASA Langley Research Center, 
Hampton VA, Jan. 1981. 

13. Kwakernaak, H., and Sivan, Linear Optimal Control 
Systems , John Wiley and Sons, New York, 1972. 

14. Safanov, M.G., and Athans, M. , "Gain and Phase Margins 
of Multi-Loop LQG Regulators ," IEEE Trans. Auto. Control , 
Vol. 22, pp 173-179, April 1977. 

15. Doyle, J.C. , and Stein, G., "Multivariable Feedback 

Design: Concepts for a Classical Modern Synthesis," 

IEEE Trans. Auto. Control, Feb. 1981. 


- 48 - 



Publications Issued During the Course of this Research 


Active Flutter Suppression Using Eigenspace* 
and Linear Quadratic Design Techniques 

William L. Garrard and Bradley S. Liebst 

University of Minnesota 
Minneapolis, Minnesota 


Abstract 


Eigenspace (ES) and Linear Quadratic (LQ) techniques are 
used to design an active flutter suppression system for the 
DAST ARW-2 flight test vehicle. The performance of the ES and 
LQ controllers are very similar in meeting control surface 
activity specifications. The ES controller provides reduced 
wing root bending moment and shear but torsional stress is 
slightly higher than with the LQ controller. The ES controller 
also results in improved flutter boundaries compared with the 
LQ controller. The LQ controller exhibits significantly better 
phase margins at the flutter condition than does the ES 
controller but the LQ design requires large feedback gains on 
actuator states while the ES does not. This results in reduced 
overall actuator gain for the LQ design. 


*This paper has been submitted to the AIAA Journal of Guidance , 
Control, and Dynamics. 
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Flexure Mode Shapes 
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(d) mode 4 , frequency « 31.7 hz, generalized 
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(j) mode 10, frequency ■ 80.3 hz, generalized mass = 0.696 lb-sec'/in 


Appendix B 

Program Descriptions and Listings 


FLUTTER 


This program generates the coefficient matrices of 
the first order state equation 
x = Ax + Bu + Tn 

This first order form is derived from the second order form 


[M xx + q a 2 ( ^> 2)x + tc s + q + [K s + q A*Jx 


+ i = l + l«xu + q + q A ? ( ^7 )G 


, - .6 , c . 2 6 . - . 6 , c . 6, - J5 „ 

+ q A 2 ( 2v ) y + q A, ( \t + 3 A n ^ 


1 2 V V 


0 V 


and i=l,2,...L aerodynamic lag states 


• 2V x — n ~ £ 

y. = -I (~) K . y . +D. X + D U + D. £ 
5 i i i xV 


or 


Mx + Cx + Kx + X, y. + Pu + Qu + Ru + SS 
i = l l 

+ T<5 + US = 0 


and 


• 2V A 

y. = -I ( — )K.y. + Dix + Eiu + Fi6 

1 c 11 


where 

M = M xx + i Aj (^) 2 

C = c s + q A x (^r) 
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K = k s + q a* 

p - % + 

Q - q A i 2v 
p = i A o 

s = i»2 <#7> 2 <7> 

T = 5 A 1 <&> 4> 

“ - 5 a o <?> 

_ . v 
Di= D. 

l 

Ei= dY 
i 

Fi= D J(i) 

x = modal coordinates and rigid body modes 
u = control surface deflections 
6 = vertical gust velocity 
V = forward velocity 
q = *5PV 2 


c = mean aerodynamic chord 

Where the actuator-aileron dynamics are 

A 

u = Ju 
u = Gu + Hu 


where 


u 



= commanded control inputs 


u = 


r — ^ 


U 


O 


U . 
1 


control surface deflections 
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And the gust model is 

2 = -(^) 2 <$ - 2(|)z - 2 . 4 6 4 ( ~) 2 n 

6 = 2 + 1.732(~)n 

where z = an intermediate gust state 
£ = characteristic gust length 
r\ = white driving noise 


Thus, if 


x = 


X 

*1 

y 2' 

y L 

u 

& 

,z . 


then neglecting SS 


and Pu 


r = 


0 

-M _1 T (1.732) (j) 
FI (1.732) (j) 

F2 (1.732) (j) 

; 

FL (1.732) (|) 

0 

(1.732) (j) 
(-2.464) (|) 2 
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I 

0 

0 . . 

0 

0 

0 

0 

-M _1 C 

-M' 1 

-M -1 . . 

. -M _1 

-M _1 (RJ+QJG) 

-M _1 U 

-M- 1 ' 

D1 

-I^ Kl 
c 1 

o • • 

0 

El(JG) 

0 

FI 

D2 

0 

2V 

K ‘ 1 
c z 

0 

E2(JG) 

0 

F2 

DL 

0 

0 • • 

2V 

• -I— K 
c L 

EL(JG) 

0 

FL 

0 

0 

0 . . 

0 

G 

0 

0 

0 

0 

0 . • 

0 

0 

0 

1.0 

0 

0 

0 . . 

0 

0 

-(V 

V 



B 


-M~ 1 QJH 
El (JH) 
E2 (JH) 

EL (JH) 
H 
0 
0 
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The open loop eigenvalues are also output. Open loop 
eigenvectors can be output as well. 

MODAL 

This program finds the optimal psuedo-inverse solution 
to the feedback gains (K) that achieve desired closed loop 
eigenvalue-eigenvector pairs (see Section on ES theory). 

The closed loop eigenvalues are also output. Closed loop 
eigenvectors can be output as well. 

COV 

This program finds the state covariance matrix for a 
given driving noise intensity. The rms aileron deflections 
and rates are determined as well. 

LOADS 

This program calculates the rms shear, bending moment, 
and torque at various station along the wing span. 

FEEDBKE 

This program calculates the closed loop eigenvalues at 
off design flight conditions. 
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O1070C V(N,H)= A<H„N) TRANSPOSED 

01030 SUBROUTINE HATTRA1 A, V j,H ? N ) 

01 090 DIMENSION A(M,N) ,0(N,M) 

01100 DO 14 1=1 ,« 














OOIOO PROGRAM LOADSUNPIJT, OUTPUT, TAPE'5= [MPIJT, TAPE/4=(1liTPIiT,COMST, 

O0110+TAPE3=C0VST,LHAT,TAPE2=LMAT) 

OO 1 2 OC* ***************************** ****t*****t***-%t*tt*f* •/**•*■ ******* 
O0110C THIS PROGRAM INPUTS CLOSED LOOP STATE COVARIANCE MATRIX 
OOt 40C (OUTPUT FROM COV), LOADS INFLUENCE MATRIX AND OUTPUTS 








2 CONTINUE 
RETURN 
END 
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